Giant van der Waals and Casimir interactions via transmission lines 
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We show that the van der Waals (vdW) and Casimir forces between neutral particles can be 
enhanced by many orders of magnitude upon changing the character of the force-mediating vac- 
uum photon-modes. By considering two induced dipoles in the vicinity of any standard electric 
transmission line, where photons can propagate in one dimension (Id), we find analytically that the 
interaction scales nontrivially with the inter-dipolar distance, resulting in much stronger and longer- 
range interaction than in free-space. This may have profound implications on the non-additivity of 
vdW and Casimir interactions in many-particle systems, and open the door for Casimir phenomena 
in Id. We discuss the possibilities of measuring this effect, e.g. in a coplanar waveguide line. 

PACS numbers: 03.70.+k, 12.20.-m, 42.50.-p, 31.30.J.- 



Introduction. Quantum electromagnetic fluctuations 
induce forces between neutral particles. This was first 
shown by London [I], to be a source of the van der 
Waals (vdW) interaction between two induced dipoles. 
Casimir and Polder [2] addressed the same interaction 
as being mediated by virtual photons from the vacuum. 
For inter-particle distances r shorter than their dipolar- 
transition wavelength A (the vdW limit), they found that 
the interaction energy falls off as 1/r 6 , whereas for longer 
distances (the retarded, Casimir limit), it scales as 1/r 7 . 
Vacuum-mediated interactions that scale differently from 
the 1/r 6 and 1/r 7 power laws are known to exist be- 
tween two polarizable bulk objects rather than point-like 
dipoles. Casimir considered the interaction between two 
parallel metal plates due to the vacuum (zero-point) en- 
ergy and found an attractive force that scales as 1/d 4 , 
with d being the distance between the plates [3]. Since 
then, many related effects have been studied, mainly 
the interactions between dielectric or metallic objects of 
various geometries [H [5]. Such experimental [6HH] and 
theoretical [T0l - fT4] studies demonstrate the role of non- 
additivity and retardation of these dipolar interactions. 

A key point in determining how these effects depend 
on distance is the space dependence of the propagation 
and scattering of the virtual photon modes that medi- 
ate the interaction. In turn, these modes depend on the 
geometrical shape of the interacting objects. Here, we 
take a somewhat different approach towards the geome- 
try dependence of vdW and Casimir-related phenomena. 
Instead of considering the interaction energy of extended 
objects with different geometries, we revisit the origi- 
nal Casimir-Poldcr configuration of a pair of point-like 
dipoles while changing the geometry of their surrounding 
environment. More specifically, we consider the energy of 
the interaction between two point-like dipoles, mediated 
by vacuum photon modes along Id. The resulting attrac- 
tive interaction is found to be much stronger and longer- 
range than its free-space counterpart, where it surpris- 



ingly scales as const. + (r/A) ln(r/A) and as 1/r 3 for short 
and long distances compared to A, respectively. This 
enhancement implies a drastic modification of Casimir- 
related effects for many-body and bulk systems in such 
a Id geometry. 

The ability to change the geometry of photon vacuum 
modes, is widely used in quantum optics, e.g. for the 
enhancement of spontaneous emission |15l 116] and reso- 
nant dipole-dipole interactions [17] . Here we consider Id 
photon modes in structures that support the propaga- 
tion of transverse-electromagnetic (TEM) modes, namely 
modes where the propagation axis, the electric and the 
magnetic fields, are perpendicular to each other. These 
are typically the fundamental transverse modes of elec- 
tric transmission lines (TLs), i.e. waveguides comprised 
of two conductors as shown in Fig. 1. They possess a 
dispersion relation w^ = \k\c and an electric field-mode 
function IT51 
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where z is the propagation axis, L its corresponding 
quantization length, A{x, y) the effective area in the 
transverse xy plane, e,j =x ^ y the polarization unit vec- 
tor, k and Wfc the wavenumber in the z direction and 
the angular frequency, respectively, and c = 1/^/JIe the 
phase velocity, fj, and e being the effective permeability 
and permittivity of the TEM mode, set by its geometry 
and materials. Unlike modes of other waveguides, such 
as optical fibers or hollow metallic waveguides, here the 
effective area A(x, y) is independent of frequency. Hence, 
considering its dispersion relation, the TEM mode forms 
an effective plane wave in Id. TLs support also higher- 
order transverse modes, such as the transverse electric 
(TE) and transverse magnetic (TM) modes in a coax- 
ial TL [18] . However, as was recently shown for metallic 
waveguides [19] , they do not contribute to the long-range 
interaction, as explained below. The existence of TEM 




FIG. 1: Geometries of transmission-line mediated vdW and 
Casimir interactions, (a) Coaxial line: two concentric metallic 
cylinders, the inner one with radius a and the outer (hollow) 
one with radius b. Two dipoles represented by black dots are 
placed in between the cylinders along the wave propagation 
direction z. They interact via modes of the coaxial line that 
are in vacuum, giving rise to a vdW-like interaction energy, 
(b) Coplanar waveguide: similar to (a). Here the central con- 
ductor of width 2a is separated from a pair of ground plane 
conductors 2b apart. 



modes as in Eq. (TTT) is correct exactly for TLs made from 
perfect conductors, however it makes an excellent ap- 
proximation of the propagation of fields in low-loss TLs 
that are in practical use [IB] . 

Analysis principles. The calculation of the interaction 
energy contributed by the TEM mode is essentially equiv- 
alent to that mediated by the vacuum in Id. We calcu- 
late it analytically by two different methods. In the first, 
we adopt the quantum electrodynamics (QED) perturba- 
tive approach [50]. We consider two identical atomic or 
molecular dipoles, with a ladder of excited levels {|n)}, 
both in their ground state \g), coupled to the vacuum 
of photon modes given by Eq. (fTl), via the interaction 
Hamiltonian 



Hj 
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Here gkjvn = v^H" d ™v ' u kj( r v) is tlie dipolar coupling 
between the dipole v — 1,2 and the kj field mode, for the 
\g) to \n v ) transition with dipole matrix element d nv , a^j 
being the lowering operator for the kj mode. The vdW 
energy is then obtained by perturbation theory as the 
fourth-order correction to the energy of the ground state 
\G) = |3i,52,0), where |0) is the vacuum of the photon 
modes 120 , 
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Here \Ij) are intermediate (virtual) states (of the free 
Hamiltonian) and E m is the energy of the state \m). The 
diagram in Fig. 2(a) illustrates one of 12 virtual processes 
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FIG. 2: Calculation methods of interaction energy, (a) QED 
perturbation theory: one of 12 possible processes (diagrams) 
that contribute to the energy correction of the state \G) = 
I <?i , <72 , 0) , Eq. pj. Here the intermediate states are |Ji) = 
|m,ff2,lfcj), I/2) = \gi,g 2 ,lk 3 lk'j'} and \I 3 ) = \gi,ria,l h >j>), 
where | lfcj) = aL|0). (b) Scattering of vacuum fluctuations: 

the vacuum field Eo(r) exists in all space and interacts with 
both dipoles at their positions ri and r 2 , hence it is also 
scattered by the dipoles. The scattered field from one dipole, 
e.g. 1, EL (r), arrives at dipole 2, resulting in an interaction 
interpreted as the vdW/Casimir interaction U. 



that contribute to this energy correction (see analysis 
details below). 

A more transparent approach towards this calcula- 
tion is that of classical electrodynamics in the pres- 
ence of an external fluctuating (quantum) vacuum field 
[3]. As shown in Fig. 2(b), we consider two point 
dipoles with polarizabilities 0:1,2(0;), subject to a fluc- 
tuating (vacuum) field Erj. The electromagnetic en- 
ergy of, say, dipole 2, is given by a sum over all k 
of -(l/2)a 2 K)[E fe (r 2 )E t fc (r 2 ) + EJ[(r 2 )E fe (r 2 )], where 
Efc(r 2 ) is the k mode component of the electromagnetic 
field at the location of this dipole. The field at r 2 includes 
two components: Eo i fc(r 2 ), the external fluctuating (vac- 
uum) field, and E^ c ' fc (r 2 ), the scattered field, which to 
lowest order is that scattered from dipole 1, where it is 
driven by the external (vacuum) fluctuations Eo,fe(ri) at 
ri. Since the TEM field exhibits diffraction-less propa- 
gation in Id, this scattered field is found by essentially 
solving 
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where al""' — a v /A(x Vl y v ) (see analysis details below). 
The vdW/Casimir energy is then obtained as the inter- 
action energy between the dipoles, which to lowest order 
becomes 



U=-2J2 
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Here we average over the quantized fields with the vac- 
uum state |0), and multiply by 4 to account for both 
field polarizations j and for the energy at dipole 1 as 
well. From Eq. Q it is clear that the scattered field 
is proportional to Eo ) fe(^i), such that the vdW/Casimir 
energy, Eq. ([5| depends on the correlation between the 
vacuum fluctuations at z\ and z<i. 

Both calculation methods yield the same result for the 
TEM-mediated interaction energy (see analysis details 
below). Here, we present the resulting energy for the 
case where one excited level |e) with energy E e and cor- 
responding wavelength A e = 2ithc/E e , has a dominant 
dipole transition, such that all other excited levels \n 7^ e) 
can be neglected, 
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where z = \z% — Z\\, is the inter-dipolar distance in the 
TL propagation direction, A v = A(x v , y u ), djr is the pro- 
jection of the dipole matrix element d e on the trans- 
verse xy plane, c.c. stands for complex-conjugate and 
Ei(x) = — P/_ dte^ t /t is the exponential integral func- 
tion. This interaction is attractive and its dependence 
over the distance z is described by the monotonously de- 
creasing function F(z), plotted in Figs. 3(a) and 3(b). 
From F(z) we can get the energy dependence at short 
(vdW regime) and long (retarded, Casimir regime) dis- 
tances, by expanding it for small and large z/\ e , respec- 
tively, 
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where 7 = 0.577 is Euler's constant. 

Predictions. The above expressions for the interac- 
tion energy in the vdW and Casimir regimes, suggest the 
possibility of a much stronger interaction than its free- 
space counterpart. In the nonretarded vdW regime, it 
decreases very slowly with z compared to the familiar 
1/z 6 scaling, whereas in the retarded Casimir regime, it 
falls off with a power-law which is four powers weaker 
than the 1/z 7 Casimir-Polder interaction. 

Let us make this comparison more quantitative by con- 
sidering a general TL with separation a between its two 
guiding conductors. It is then reasonable to assume that 
the effective TEM mode area A(x,y) scales as a 2 , a be- 
ing the relevant length-scale. For example, in a coax- 
ial line [see Fig. 1(a)], by normalizing the electric field 
of the TEM mode, we find ^/A = ^2Trln(b/a)p, with 
p = y/x 2 + y 2 [IS]. Since a < p < b, taking e.g. 
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FIG. 3: The TEM-mode mediated interaction potential as 
a function of inter-dipolar distance z. (a) Log-log plot of 
F(z), Eq. (pSl. For long distances, z ^> A e , the linear de- 
pendence implies a power-law behavior as in Eq. (7b. (b) 
F(z) at short distances, (c) Log-log plot of the ratio between 
TEM-mediated energy U(z), Eq. (fSJ) , and the free-space vdW 
energy in short distances z < A e , U fs {z) — - 4g -/ e 2 fie ^r [20] 
, with jd^-| 2 = (2/3) I d e | 2 and ^J A 1>2 = a (see text). Here 
a = 10~ 4 A e is taken, consistent with typical cases consid- 
ered later, where a ~few /im and E e /(2irh) ~few GHz. From 
z ~ 10~ 3 A e and beyond, the huge enhancement of the interac- 
tion w.r.t its free-space counterpart is apparent, (d) Same as 
(c), but for long distances z>A e where the free space energy 

takes the Casimir-Polder form, U fs (z) = ~ ' d ' ' 
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b = 2a and p = a we obtain V A ~ 2.1a. Moreover, 
the energy in Eq. (|6| depends on the projection of the 
dipole element on the transverse plane. Assuming the 
polarizability is isotropic on average, we can estimate 
|d^-| 2 = (2/3)|d e | 2 (Supplementary Information). Using 
this estimation with WAi^ = a, and taking the reason- 
able values of a/X e < 1, we calculate the ratio between 
the TEM-mediated energy and its free-space counterpart. 
This ratio is plotted in Figs. 3(c) and 3(d) for short and 
long distances respectively, using the known free-space 
expressions J2HJ- From distances as short as z — 10~ 3 A e , 
the TEM-mediated interaction becomes orders of magni- 
tude stronger, where this ratio increases drastically with 
z. 

At very short distances, however, the free-space vdW 
interaction 1/z 6 is stronger. This occurs at z < a, where 
the dipoles are close enough such that they do not " sense" 
the TL structure. The transition to free-space behavior 
in this regime can be described by including higher-order 
transverse modes. These modes become significant at 
such short distances, where their contributions sum up 
to give the free-space result. In the Supplementary In- 
formation we calculate, for a coaxial TL, the interaction 
energy due to the TE and TM modes. We find that the 



energy contribution of the TEi m and TMi m modes with 
cutoff frequency ck\ m scales like K (ki m z) and e~ klmZ re- 
spectively, where Kq{x) is the modified Bessel function. 
Since fc; m > n/a, at long distances z > a both decay ex- 
ponentially and are negligible with respect to the TEM 
mode energy, however at short distances z < a, we nu- 
merically verify that the dominant TM modes sum up to 
give exactly the free-space interaction (Supplementary 
Information). This was also recently shown for the dis- 
persion interaction in a metallic waveguide [19) . and we 
expect it to hold for all TLs. Namely, for distances z ^> a 
we can indeed consider only the TEM mode, whereas the 
free-space result is recovered for z <C a, owing to the 
other transverse modes. 

It is important to consider how the predicted effects 
can be measured. Here we focus on the coplanar waveg- 
uide (CPW) TL gl] [Fig 1(b)] that is extensively used 
in the emerging field of circuit-QED 22-26J, and in 
which the dynamical Casimir effect has been recently 
demonstrated [27]. We consider the dominant dipole 
transition of a pair of superconducting transmons [24] 
at distance z, both capacitively coupled to the CPW. 
Then, the TEM-mediated interaction should induce a z- 
dependent energy shift on the dipole-transition levels, 
U(z) from Eq. (pi). In order to estimate this energy 
shift, we refer to a recent experiment [25], where the 
dipole frequency is E e /{2 / Kh) ~ 5GHz and the dipole 
coupling to a closed CPW cavity of length A e , can reach 

g ~ it x 720MHz. Using g = \/^P^^= and L = A e 

[251 l2"o] . we can extract Id^l/v^ for each dipole, ob- 
taining [U(z)/F(z)]/h ~ 0.84MHz, where h is Planck's 
constant. E.g., for distances z = 0.001A e or 0.01A e (both 
much larger than a ~few /im), the energy shift becomes 
U(z)/h sw 1.8 or 2.47MHz respectively, about twice the 
dephasing rate of the dipole, 1MHz [25] . that limits the 
resolution of the shift. However this resolution can be 
considerably improved, as in Ref. [26], where a dephas- 
ing time of about 20/is is reported. Then, by taking 
E e /h ~ 2GHz with the parameters of Ref. [23], one 
can obtain U{z)/h s=s 28MHz for z = 0.01A = 1.5mm, 
which is much larger than the dephasing. Probing the 
interaction in the retarded regime, where z > A e , is 
currently more challenging: for the latter case, with 
z = 2A e = 30cm, we obtain U(z)/h w 6.62KHz, which is 
too small to be observed, although future improvement 
may make this regime observable as well. Nevertheless, 
these results imply the remarkable possibility to directly 
observe the vdW and Casimir interaction between a sin- 
gle pair of point-like dipoles (here the size of the dipoles 
is about l/im) over a wide range of macroscopic distances 
where the interaction scales non-trivially as in Eq. d6]). 

Moreover, since the open geometry of the CPW al- 
lows for the coupling of the TEM mode to clouds of 
trapped Rydberg atoms or polar molecules above its sur- 
face [25] [25], the possibility exists to explore the predicted 



long-range interaction in a many-body setting. This may 
entail a modification of the non-additivity of the vdW 
and Casimir interactions, which currently attracts con- 
siderable interest [5J. In free-space, the vdW energy of a 
gas is approximately additive, namely it can be obtained 
by pairwise summation of the vdW energies of all pairs, 
as long as a/r 3 is small, with r a typical inter-dipolar sep- 
aration 4]. This scaling is however a direct consequence 
of dipole-dipole interactions in free space, and is expected 
to change when the, e.g. atomic dipoles, interact also via 
the TEM mode of the CPW, for the same reasons that 
the 1/r 6 scaling was shown here to change. Such atomic 
clouds are anticipated to exhibit non-additive effects at 
smaller densities than usual, which in turn may influence 
their thermodynamic properties as well as their effective 
dielectric constant. Furthermore, it will be interesting 
to consider how the long-range interaction we predict is 
modified at finite temperatures. 

Analysis details. Beginning with the QED perturba- 
tive approach, the sum from Eq. (13j) includes 12 different 
terms, each corresponding to a different set of interme- 
diate states and represented by a diagram [2UJ. E.g. the 
term corresponding to the diagram from Fig. 2(a) is 
given by 
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Summing all 12 terms and then performing the integra- 
tion over kl we arrive at 
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(9) 
where k n — E n /(hc). Taking a single excited level |e) 
for each dipole, the integration on k is performed (using 
regularization) and Eq. ^ is obtained. 

Turning to the calculation in the vacuum scattering ap- 
proach, we observe that the scattered field is proportional 
to the Green's function of the Helmholtz equation in Id, 
Eq. (Ik, which is found to be 2lk\e l ^ z ~ Zl ^ Then, insert- 
ing this field into the energy Eq. ([5| , where the quantum 

(vacuum) field in Id is Eo,jy(z) = i\l ^r-jy e lkz ej<ikj 
[the index j is suppressed in the Eqs. Q, (pM)], we arrive 
at 



U = 



lie 



2ire 2 A 1 A 2 



dka 1 (k)a 2 (k)k 2 sm(2kz). (10) 



By assuming the polarizabilities ai,2 of a system with 
a discrete set of transitions as in an atom, a(k) = 

I J2 n E2-h 2 c 2 k 2 ESI' anci tnen transferring the integra- 
tion to the imaginary axis in a complex fc-plane (Wick 



rotation, with regularization), we arrive at a result equiv- 
alent to that in Eq. ^: it is described by a Meijer G- 
function, that allows to analytically obtain exactly the 
same results as in Eq. Q when |d^-| 2 = (2/3)|d e | 2 is 
assumed, and can be numerically calculated for a general 
z. 

Conclusions. To conclude, we have studied how 
transmission-line structures, that are typically used to 
transmit electromagnetic signals in electronic devices, 
can in fact transmit vacuum fluctuations between dipoles 
and drastically enhance dispersion forces. The key point 
is that the modes of the electromagnetic vacuum are 
changed in the presence of TLs and include the TEM 
mode that propagates in Id. We have shown how the re- 
sulting vdW and Casimir-Polder interactions can become 
longer- range and larger by orders of magnitude than their 
free-space counterparts. To this end, we have analyt- 
ically found, via two independent methods, an expres- 
sion for the dominant TEM-mediated interaction at all 
inter-dipolar distances, Eq. (pi), and described how the 
free-space result is restored at very short distances, by 
including higher order modes in our calculations. Our 
approach assumed that the TL is comprised of perfect 
conductors at zero temperature. This is justified for a 
coplanar-waveguide line used in low-temperature, super- 
conducting circuit-QED experiments, for which we have 
estimated that the effect may be directly measured. We 
stress that unlike other studies, where the circuit-QED 
system is used as a single-mode cavity [211 EH [30] , here 
we consider the open waveguide geometry, which allows 
for quantum optics in Id |23| . 

Our result may pave the way towards further intrigu- 
ing Id Casimir effects upon considering extended inter- 
acting objects and many-particle systems that are cou- 
pled to a transmission line. These configurations may re- 
veal more complex Casimir phenomena than the simple 
dipole-dipole non-retarded vdW interaction due to two 
major effects; i.e., retardation and non-additivity [S]. As 
discussed above, both of these may be drastically modi- 
fied in a Id geometry, namely, by the presence of a trans- 
mission line. 
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